GLEARiMGKOUSE  FOR  FEDERAL  SCIENTIFUC  AND  TtCKNSGAL  mfi  iMATiON  CFSTI 

DOCUMENT  MANAGEMENT  BRANCFl  410.1  i 


LIMITATIONS  IN  KEPRODUCTiQN  QUALITY 


ACCESSION  If  41^  4ovyy7 


n  I .  llfE  REGRET  THAT  LEGIBILITY  OF  THIS  DOCUMENT  IS  IN  PART 

UNSATISFACTORY.  ncPRCDUCTIQN  HAS  BEEN  MADE  FROM  BEST 
AVAILABLE  COPY. 


2.  A  PORTION  9F  THE  ORIGINAL  DOCUMENT  CONTAINS  FINE  DETAIL 
WHICH  MAY  MAKE  READING  OF  PHOTOCOPY  DIFFICULT. 


n  3.  THE  ORIGIKAL  DOCUMENT  CONTAINS  COLOR.  BUT  DISTRIBUTION 
COPIES  ARE  AVAIUBLE  IN  BLACK~AND>WHh  £  REPRODUCTION 
ONLY. 


n  4.  THE  INITIAL  DISTRIBUTION  COPIES  CONTAIN  COLOR  WHICH  WILL 
BE  SHOWN  IN  BLACK-ANO-WHITE  WHEN  IT  IS  NECESSARY  TO 
REPRINT. 


n  S.  LIMITED  SUPPLY  ON  HAND:  WHEH  EXHAUSTED.  DOCUMENT  WILL 
BE  AVAILABLE  IN  MICROFICHE  ONLY. 


n  I.  LIMITED  SUPPLY  ON  HAND:  WHEN  EXHAUSTED  DOCUMENT  WILL 
NOT  BE  AVAILABLE. 


Q  2-  DOCUMENT  is  AVAILABLE  IN  MICNOFICHE  ONLY. 

Q  8.  DOCUMENT  AVAILABLE  ON  LOAN  FROM  CFSTI  (  TT  DOCUMENTS  ONLY). 

□  » 


PROCESSOR: 

0 


m  B/S4 


BEST 

AVAILABLE  COPY 


604447 


i»F:RL**64-608 


?!?E  REFRACTED  ARRIVAL  FROM  A  LAYER 
T.  W.  Spencer 

California  Research  Corporation 
Jj&  Habra,  California 


ContrAct  AP19 (60^1 )  ”8^4^ 

« 


SclenMflc  Report  Number  6 


June  15 »  196^ 
Project  8652 
Task  865205 


7^;: 


HARO  copy 

microfiche 


Prepared  for 


AIR  FORCE  CAMBRIDGE  RESEARCH  LABORATORIES 
OFFICE  OF  AEROSPACE  RESEARCH 
UNITED  AIR  FORCE 

BEDFORD,  MASSACHUSETTS 


WORK  SPONSORED  BY  ADVANCED  RESEARCH  PROJECTS  AGENCY 

PROJECT  VELA-U]«IFORM 

ARPA  Order  No.  I8O-6I 
Project  Code  No.  8IOO  Task  2 


OOCIRA  C 


D  DC 

rhJ2i2EJILflC^ 
StPl  1964  I 

OfegglLU  \I  IJbiE 


Requeot  for  additional  copies  by  Agencic-rc  ui  the 
Department  of  Defense,  their  contractors,  and  other  gdvemraent 
agencies  should  be  directed  to  the: 

DEFENSE  DOCUMENTATION  CENTER  (DDC) 

CAMERON  STATION 
ALFJaNDRIA,  VIRGINIA 

Department  of  Defense  Contractors  rau-st  be  established  for  DDC 
services  or  have  their  ''need-to-kno'-f”  certified  by  the  cognizant 
military  agency  of  their  project  or  contract. 

All  other  persona  and  organizations  should  app^’y  to 

the: 

U.  3.  DEFARTMENr  OF  COMMERCE 
OFFICE  OP  TECHNICAL  SERVICES 
WASHINGTON,  D.  C.  20230 


ill 


THE  REFRACTED  ARRIVAL  FROM  A  LAYER 
T.  W.  Spencer 

ABSTRACT 

Hlgh-freqiiency  geometric  ray  theory  is  used  to  Investi¬ 
gate  the  refracted  arrival  from  a  high-speed  layer  embedded  in  an 
infinite  medium.  The  effect  of  changing  the  layer  thickness  to 
dominant  wavelength  ratio  (E/  )  and  the  range  to  depth  ratio 

is  analyzed  for  a  point  corapresslonal  source.  The  results 
approximate  the  exact  solution  when  E/)^sE:  !  .  The  theory  predicts 
shingling  and  shoots  that  it  is  distance-limited.  Tlie  maximum 
distance  at  which  shingling  is  observed  decreases  as 
decreases.  Beyond  a  certain  minimum  distance,  the  refracted 
arrival  consists  of  a  superposition  of  waves  whose  phase  veloci¬ 
ties  appi-oach  the  compressional  velocity  in  the  layer  at  large 
distances.  These  waves  include  the  reflected  waves  which  cross 
the  layer  at  least  once  in  the  compressional  mode  and  the  head 
waves  which  are  generated  by  reflected  waves  which  do  not  cross 
the  layer  in  the  compressional  mode.  There  are 

reflected  ray  paths  which  cross  the  layer  h  times  in  the  com¬ 
pressional  mode  andrn  times  in  the  shear  mode  vdilch  arrive  at 
exactly  the  same  time.  A  formula  is  derived  which  expresses  the 
composite  amplitude  of  these  rays  in  the  fonn  of  a  series  which 
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contains  either  m  ur  fl+l  terras  (whichever  is  smaller).  The 
travel -time  ''urves  for  all  I’ays  with  the  same  !T1  value  approach 
the  same  asymptote  regardless  of  the  H  values.  At  large 
distances  the  travel-time  curves  for  waves  whose  fTl  values  differ 
by  one  approach  asymptotes  which  differ  by  .  At  large  dis¬ 
tances  and  for  thin  layers,  the  total  amplitude  obtained  by 
sumntlng  all  rays  with  the  same  HI  value  satisfies 

The  depend  only  on  the  density  and  velocity  ratios.  Because 
the  head  waves  decay  like  (P/H)  ,  the  reflected  waves  predomi¬ 

nate  at  large  distances.  If  the  duration  of  the  direct  ware  is 
less  than  ,  the  refracted  arrival  consists  of  a  series  of 
events  -  each  associated  with  a  different  fD  value.  At  large 
distances,  the  onsets  of  the  events  are  separated  by  ,  the 
amplitudes  decay  like  (^/H)  and  tlie  ’'iffj  determine  the  relative 
amplitudes . 


V 


CCNTEI'ITS 

Page 


INTRODUCTION  1 
THE  EMBEDDED  LAYER  3 
DEGENERACY  4 
REMOVAL  OF  DEGENERACY  6 
HIGH-FREQUENCY  THEORY  10 
CHARACTERISTICS  OF  THE  MUI.TIFLY  REFLECTED  WAVES  19 
THE  REFRACTED  ARRIVAL  25 
CONCLUSIONS  33 
ACKNOWLEDGMENTS  37 
REFERENCES  38 
FIGURES  40 


INTRODUCTION' 


The  refracted  wave  generated  at  the  plane  interface 
between  two  half  spaces  has  been  studied  extensively  by  Heelan 
(1953)^  Strick  (1959),  Dix  (1961)  and  Gilbert  and  Laster  (1962). 
For  studies  of  explosion-generated  waves,  the  results  for  a 
point  corapressional  source  are  of  particular  significance. 
High-frequency  theory  predicts  that  the  refracted  varies 

with  time  in  a  manner  which  is  determined  by  the  integral  of  the 
direct  wave  (measured  at  distances  larj?;e  compared  with  the  wave¬ 
lengths  of  interest)  and  the  amplitude  decays  with  distance  like 

x-l/2  Jill 

y  ^  is  the  horizontal  source-receiver  separation  and 

L  is  the  distance  traveled  in  the  refractor.  Donato  (1964a) 
has  determined  the  minimum  range  (in  terms  of  the  dominant  wave¬ 
length)  at  which  the  high-frequency  approximation  is  valid. 

When  the  high-velocity  medium  Is  finite  in  thickness, 
the  problem  may  be  considerably  more  complicated.  To  demonstrate 
this,  let  T  be  the  duration  of  the  head  wave  for  an  infinitely 
thick  layer  and  let  be  the  head  wave  travel  time.  The  results 
for  an  infinitely  thick  layer  are  not  valid  if  waves  which  are 
multiply  reflected  within  the  layer  or  other  wave  types  arrive 
in  the  interval  ^  "tj..  +  T.  Multiple  reflections  distort 

the  refracted  arrival  when  T  is  large,  the  layer  is  thin  or  the 
offset  is  large.  When  any  of  these  conditions  are  satisfied, 
the  individual  phases  are  not  resolved  and  the  character  and 
amplitude  of  the  refracted  arrival  are  determined  by  the  way  in 
which  the  individual  phases  Interfere  with  one  another.  Inter¬ 
ference  explains  the  amplitude  dependence  on  offset  and  on  the 
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dominant  wavelength  to  thickness  ratio,  the  discrepancy  between 
the  actual  velocity  in  the  layer  and  the  measured  velocity,  the 
error  in  determining  the  tr-ue  beginning  of  the  arrival  and  the 
phenomenon  of  shingling^. 

Dunkin  (1963)  and  Donato  (l^j-o^b)  have  investigated 
theoretically  the  first  extremum  produced  by  the  interference 
of  the  head  wave  and  the  compresslonal  wave  reflected  off  the 
bottom  of  the  layer.  Their  theoretical  predictions  for  the 
dependence  of  phase  velocity  on  layer  thickness  and  the  depend- 
once  of  amplitude  on  range  are  in  general  agreement  with  model 
experiments  (Lavergne,  1961,  and  Levin  and  Ingram,  1962). 
Various  lines  of  experimental  evidence  (e.g.,  shingling)  indi¬ 
cate  that  in  certain  situations  the  first  extremum  may  be  so 
weak  thac  it  is  not  detected.  In  this  paper  high-frequency 
geometric  ray  theory  is  used  to  compute  the  complete  refracted 
wave  train. 


■^Shingling  occurs  when  the  phase  ocity  exceeds  the  group 
velocity.  As  the  range  Increases,  peaks  and  troughs  move 
forward  through  the  envelope  which  defines  the  refracted 
arrival.  In  this  process,  the  amplitude  of  the  first  extremum 
decreas(js  and  it  is  eventually  lost  in  the  noise.  At  this  off¬ 
set,  a  later  extremum  is  selected  to  define  the  time -distance 
curve.  At  each  offset  where  an  extremum  is  lost,  there  is  a 
discontinuity  in  the  time-distance  curve  and  a  new  shingle  is 
added  corresponding  to  a  later,  larger  amplitude  extremum. 


THE  EMBEDDED  LAYER 


Consider  the  problem  of  determining  the  refracted 
arrival  from  a  hlgh«vel>Dcity  layer  embedded  in  an  infinite 
medium.  This  model  avoids  the  complications  Introduced  by 
putting  in  a  free  surface.  The  significance  of  the  geometric 
quantities  is  indicated  in  Figure  1.  The  top  of  the  layer  is 
at  8  depth  H  below  the  soui?ce  and  at  a  depth  Z  below  the 
receiver.  The  radiation  field  is  axially  symmetric  about  the 
7,  axis.  The  Z  axis  passes  through  the  source  and  is  normal  to 
the  Interface.  Tlie  plane  Z  =  0  coincides  with  the  top  of  the 
layer.  Our  computations  are  made  for  a  horizontal  refractor 
so  that  Z»H.  A  dipping  I'efractor  can  also  be  investigated  by 
using 

(z-HI 

to  relate  Z  and  H.  Here  the  refractor  dips  at  an  angle  (j)  with 
respect  to  the  horizontal  plane  in  which  the  source  and  receiver 
are  located  and  is  the  source-receiver  separation  in  the 
horizontal  plane.  E  is  the  layer  thickness.  V  is  the  corapres- 
slonal  velocity,  or  is  the  shear  velocity  and  b  is  the  density. 
Subscript  1  refers  to  the  infinite  medium  and  subscript  2  refers 
to  the  layer. 

The  critical  angle  for  P, (incident  compressional, 
transmitted  compressional)  is  9^  and  the  crli-ical  angle  for  RS^ 
(incident  compressional,  transmitted  shear)  is  9^  . 
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DEGENERACY 

When  we  attempt  to  take  into  account  all  the  waves 
which  are  multiply  reflected  within  the  laye  r**  I  vrc  ^uuuCvA  lately 
encounter  the  problem  of  aegeneracy.  The  simplest  case  of 
degeneracy  arises  when  the  ray  crosses  the  layer  once  in  the 
compressional  (P)  mode  and  once  in  the  shear  (r ;  mode.  In 
Figure  1  we  see  that  there  are  two  such  rays  which  arrive  at 
exactly  the  same  time  and  are  therefore  degenerate.  These  two 
rays  differ  in  the  sequence  in  which  the  P  and  S  legs  of  the 
path  are  traversed.  This  type  of  degeneracy  cannot  be  removed 
by  changing  the  layer  thickness  or  velocity.  It  is  a  conse¬ 
quence  of  the  fact  that  two  different  types  of  waves  can  propa¬ 
gate  in  an  elastic  medium.  Degene rnoy  occurs  only  for  rays 
which  cross  the  layer  in  both  modes. 

To  reduce  this  idea  to  quantitative  terms,  we  introduce 
the  following  definitions; 

r\  number  of  crossings  in  the  P  mode, 

Dfl  number  of  crossings  in  the  S  mode. 

Specifying  n  and  ITI  determines  a  family  of  degenerate  rays.  The 
number  of  rays  in  the  family  (i.e.,  the  degeneracy)  is  given  by 


(1) 
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^Even  when  n  -  6  and  m  =  2,  the  degene rscy  is  28.  If  highly 
degenerate  events  are  truly  significant,  a  way  must  be  found 
for  removing  the  degeneracy  -  otherwise  the  ray  theoi'y  becomes 
extremely  inefficient  from  a  computatlonPxl  standpoint.  Also, 
because  all  the  degenerate  rays  arrive  at  the  same  time,  only 
the  composite  event  has  jiiyslcal  significance.  For  the  embedded 
layer,  a  simple  technique  ejclsts  for  removing  the  degeneracy. 
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REMOVAL  OF  DEGENERACY 

Per  Illustrative  purposes,  we  consider  waves  which 
leave  the  source  and  arrive  at  the  detector  in  the  compressional 
mode.  The  total  amplitude  of  a  degenerate  event  must  be 
expressible  in  the  forrii 


The  T's  are  transmission  coefficients.  The  etter  subscripts 
indicate  the  mode  of  propagation  before  f..  d  after  transmission. 
The  numbers  Indicate  the  media  in  which  the  wave  pi'opagal  J. 

gives  the  total  amplitude  of  all  waves  which  make  the  first 

nrn 

and  last  crossing  In  the  layer  in  the  compressional  mode.  The 
significance  of  each  quantity  is  indicated  below: 


Symbol 

^nm 

^nro. 


First  Crossing 
P 


P 

S 


Last  Crossing 
P 
P 
S 
S 


The  transmission  coefficients  appear  explicitly  in  (2).  , 

functions  of  the  reflection  coefficients 
and  R^  .  There  are  only  four  reflection  coefficients 
because  the  wave  can  approach  the  Interfaces  only  from  within  the 
layer  and  the  external  media  have  identical  properties. 


^PP ' 


Formulas  wlxlch  give  and  in  terms 

of  the  reflection  coefficients  can  be  derived  by  using  the  pro¬ 
cedure  outlined  below  for  .  Partition  the  ni  shear  crossings 
into  groups.  Sarh  group  must  contain  at  least  one  shear 
crossing.  ^  satisfies 

1  minimum  (n- 1  ^m)  ^  ro^\^  n^2  . 

The  number  of  ways  of  distributing  m  identical  things  among 
groups  is 


/S 

If  there  are  ^ 

for  the  first  and  last  leg  must  be  in  the  P  mode).  The 

number  of  ways  of  distributing  the  n  compresslonal  crossings 
Into  ^+1  groups  is  •  Consequently,  the  total  nijunber 

of  P-S  sequences  in  which  the  S  are  arranged  in  ^  gi-oups  is 

The  exponents  for  the  reflection  coefficients  can  be 
derived  Just  as  easily.  There  are  ^  groups  of  S.  This  Indicates 
that  there  are  ^  conversions  from  P  to  S  and  from  S  back  to  P, 

Rp3  and  Rsp  enter  in  the  form 

Let  be  the  nujiiber  of  shear  Ci'osslngs  in  the  l^h  group. 

is  the  number  of  SS  reflections  in  the  1^^^  group.  The  total  number 

of  SS  reflections  is 

(%0-m-f  • 
l-j 


groups  of  3,  there  must  be  j.4-1  groups  of  P  (since 
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i3  the  exponent 
is  given  by  H  r  fO-l. 


of  Rqq .  The  total  number  of  reflections 


We  have  already  accounted  for 


+  reflections.  Therefore,  the  number  cf  PP  reflections  is 

The  final  expression  for  is 


r\sr,2  . 


The  final  expression  for  the  total  amplitude  is 


The  total  amplitude  for  the  trivial  cases  where  either  f)  or  ITl 
is  zero  is 


Instead  of  computing  the  amplitudes  of  Individual 

degenerate  rays,  we  need  evaluate  only  171  or  +1  (whichever  is 
smaller)  terms  in  a  series.  It  is  this  fact  which  permits  the 
rapid  computation  of  highly  degenerate  events.  Furthermore,  the 
use  of  (3)  avoids  the  complicated  bookkeeping  required  to  insure 
that  all  members  cf  the  degenerate  family  have  been  included. 
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In  the  derivation  we  confined  our  attention  to  a  geometrical 
situation  In  which  the  source  and  receiver  were  on  the  same  side 
of  the  layer  and  to  rays  which  leave  the  source  and  arrive  at  the 
detector  In  the  P  mode.  The  same  method  yields  the  total  ampli¬ 
tude  when  the  waves  leave  the  source  and  arrive  at  the  detector  In 
either  the  P  or  S  mode  and 

(a)  the  source  Is  in  the  half-space  and  the  receiver  Is 
in  the  layer  (or  vice  versa), 

(b)  the  source  and  receiver  are  on  opposite  sides  of 

the  layer, 

(c)  the  oource  and  receiver  are  in  the  layer. 
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HIGH-FREQUENCY  THEORY 


Knowing  we  can  derive  Information  about  the 

i:'efracted  arrival  by  using  the  method  of  generalized  reflection 
and  transmission  coefficients  (Spencer,  i960).  To  use  this 
method,  we  require  an  integral  representation  for  the  Laplace 
transform  of  the  source  radiation  field.  The  vertical  component 
of  particle  velocity  in  a  spherically  synunetrlc  compression  .1 
wave  is 


(6) 


The  source  function,  )[  ,  is  time-limited  and  merely  expands  or 
contracts  with  "tj  but  does  not  change  shape.  The  amplitude 
spectrum  of  is  assumed  to  peak  at  a  dominant  frequency 
■fj  -  i/ij .  f'  is  the  distance  between  source  and  receiver. 

At  dlstanceE'.  w.hich  are  large  compared  with  the  dominant  wave¬ 
length,  the  direct  wave  propagates  without  change  of  shape  and 
the  time  variation  is  given  by  )/ { T  ) .  CL  is  the  rt'dlus  of  the 
sji^erlcal  cavity  and  is  the  peak  pressure  applied  to  the 
cavity  wall. 

The  Laplace  transfonn  of  (5)  is 
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Using  the  method  of  generalized  coefficients  gives  the  Laplace 
transform  of  the  composite  event  specified  in  equation  (3)  as 


d>  ^ 


where 


Cagnlard’s  method  (1961)  gives  the  response  function 
corresponding  to  in  terms  of  a  contour  lutegrai.  The 
numerical  integration  is  time  consuming  even  on  a  c '>mputer  in 
the  IBM  709^  class.  For  a  line  source,  the  Bessel  function  in 
(6)  is  replaced  by  the  exponential  (  the  line  source 

coincides  with  the  7.  axis).  Garvin  (1956),  by  applying  Cagniard's 
method  to  the  line -source  problem,  obtained  a  closed  algebraic 
expression  for  the  response.  This  function  can  be  evaluated 
rapidly  on  the  IBM  709^*  The  total  solution  for  our  problem  is 
built  up  by  superposing  the  response  functions  for  all  rays  which 
arrive  during  the  time  Interval  of  Interest.  If  there  are  a  great 
many  rays,  the  computation  of  the  total  response  is  time  consiuning 
even  though  the  computation  time  per  ray  is  nominal.  Because  of 
this  and  other  difficulties,  we  have  derived  an  approximation  to 
the  exact  response  function  which  is  valid  in  some  high-frequency 
limit.  The  low-frequency  cutoff  will  become  apparent  later  when 
we  plot  the  amplitude  of  the  refracted  arrival  against  the  dominant 
wavelength  to  thickness  ratio. 
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The  high-frequency  asymptotic  expansion  of  (6}  can 
be  obtained  by  using  the  method  of  steepest  descent  and  a  modi¬ 
fication  of  the  method  of  stationary  phase  discussed  by  Erdelyi 
(1956).  'fhe  asymptotic  representation  for  the  composite  reflec¬ 
tion  is  designated  by  a  superscript  P  and  is 


=r''  CL  -jk  £  1 

{,0,1  p.  nm  ) 


where  <=f  refers  to  the  imaginary  part  of  and  to 

the  real  p£irt.  Arons  and  Yennie  (1950)  derived  a  similar  result 
by  considering  the  reflection  of  a  plane  wave  from  a  llquld-pclld 
interface.  Equation  (7)  shows  that  their  result  remains  valid 
for  a  multiply  reflected  wave  which  is  Initially  spherical  in 
shape.  is  the  angle  of  Incidence  at  the  upper  interface  of 

the  rays  which  cross  the  layer  n  times  in  the  compresslonal  mode 
and  ni  times  in  the  shear  mode.  A^^  (-lSLn@nm  )  ^  function  of 

the  plane  wave  reflection  and  transmission  coefficients.  Inside 
the  critical  angle  (  0^,  )*  Apyp(-Csm6p^  )  is  pure  real, 

the  first  term  in  (7)  vanishes  and  the  time  variation  in  the 
reflected  wave  is  determined  by  the  source  function.  Each  reflected 
wave  propagates  with  a  phase  velocity  which  depends  on  the  range. 
Waves  which  cross  the  layer  at  least  once  in  the  compresslonal 
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mode  (H— 1)  propagate  with  phase  velocities  which  approach  the 
compressional  velocity  in  the  layer  at  large  distances.  All  other 
•reflected  waves  have  phase  velocities  which  approach  either  the 
compressional  velocity  in  the  half  space  or  the  shear  velocity 
in  the  layer.  Therefore,  beyond  a  certain  minimum  distance 
(which  depends  on  the  time  Interval  of  Interest)  only  reflected 
waves  which  sat-isfj-  the  condition  Ti-ii^l  can  contribute  to  the 
refracted  arrival.  For  these  waves  the  angle  of  Incidence  is 
less  than  the  critical  angle  and  the  wave  shape  is  deteirolned 
by  tha  source  function. 

takes  into  account  the  diminution  in  amplitude 
produced  by  geometric  spreading 


is  an  Implicit  function  of  the  geometric  parameters 
through  the  relation 

±  Sine  '•nf  .4 _ ■  "'■I 


and  is  the  norvalized  travel  tline 


w 
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These  relations  show  that  except  for  the  factor  (H/Cl)  bj  V, /P^  , 
the  resp  .nse  expressed  as  a  function  of  depends  on  the 

three  geometric  parameters  p/H,  E/H  and  Z/H  and  on  the  dominant 
wavelength  to  thickness  ratio  through  the  parameter 

'bi  ll 

W^r  H  H 

The  reflected  waves  which  do  not  cross  the  layer  in 
the  compresslonal  mode  ( H  =0  )  generate  head  waves.  The  head 
waves  which  propagate  with  a  pliase  velocity  are  compo.ients 
of  the  refracted  arrival.  The  wave  reflected  from  the  top  of 
the  layer  generates  the  ordinary  head  wave  indicated  by  the 
segment  ED  in  the  upper  diagram  in  Figure  2.  This  wave  travels 
in  th'=»  'direction  normal  to  ED  with  a  velocity  Vj  and  has  a 
constant  phase  velocity,  V2  .  The  incident  wave  produces  both 
a  transmitted  compresslonal  wave  (P^P^)  and  a  transmitted  shear 
wave  (P^Sg)  in  the  layer.  These  two  wavefronts  intersect  at  the 
interface  as  long  as  the  phase  velocity  in  the  incident  wave 
exceeds  .  When  the  phase  velocity  in  the  incident  wave  drops 
below  the  wavefronts  separate.  In  this  process  an  internal 

head  wave  is  generated  (BD),  This  wave  travels  in  the  direction 
normal  to  BD  with  a  velocity  . 

The  lower  diagram  shows  what  happens  when  the  transmitted 
shear  wave  and  the  head  wave  reflect  off  the  bottom  of  the  layer. 

BQ  is  the  reflected  part  of  the  head  wave.  The  incident  shear 
wave  (PC)  generates  a  reflected  compresslonal  wave  (QP)  and  a 
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reflected  shear  wave  (PBA).  At  a  point  on  the  lower  Interface 
where  the  phase  velocity  of  the  incident  shear  wave  drops  below 
>  these  two  wavefronts  begin  to  separate  A  new  head  wave 
is  generated  whose  wavefront  coincides  with  BQ.  Each  time  the 
unconverted  shear  wave  ( h  =  0 )  is  reflected  back  into  the  layer, 
a  new  head  wave  is  generated  whose  wavefront  coincides  with  the 
wavefront  of  head  waves  generated  by  earlier  shear  reflections. 

The  high-frequency  representation  for  head  waves  which 
propagate  with  a  p^iase  velocity  is 


CJLf _  r;^  (±  .X I  ,1.) 

I-olJ  r,  H  ) 

M-m 


MWi. 


where 


(8) 


and  Pjij  and  are  defined  subsequently.  The  time  variation 

is  determined  by  the  integral  of  the  source  function.  At  large 

_2 

offsets  the  amplitude  decays  like  (  p/li)  .  The  wavefront  of  the 
head  wave  is  tangent  to  the  reflected  wave  which  generates  it  at 


the  critical  distance 


fsi, 

H 


V» 


X 


i+1 


E 


Hr  (£-Sri 


and  the  no  mall  zed  travel  time  is 
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The  head  waves  are  separated  by  equal  time  Intervals 

H/V,  H/V,  i^/  ‘ 

which  do  not  depend  on  the  offset.  This  fact  suggests  that  if 
the  refracted  arrival  consisted  of  head  waves  alone,  the  condi¬ 
tion  =  2T^  would  lead  to  constructive  Interference  when  the 
all  have  the  same  polarity. 

The  /^depend  on  the  velocity  and  density  ratios  and 
are  defined  by 


where 

is  the  complex  conjugate  of -p(f).  ^  may  be  identified 
with  the  sine  of  the  angle  of  incidence.  The  form  of -f  ) 


(10) 


depends  on  the  value  of  IT)  as  follows: 


fg)-  R^j,(0  ,  m-o  j 

. 

■f(S)  contains  square  roots  of  the  form  (.F^- -ff  ) .  A  one-to-one 
correspondence  between  points  in  the  f  and  f  (f  )  planes  is  estab¬ 
lished  by  Introducing  cuts  of  finite  length  along  the  rea]  ^  -axis 
between  each  pair  of  branch  points.  Each  square  root  is  defined 


to  be  positve  real  on  the  real  f  -axis  to  t)ie  right  of  the  cut , 
At  large  offsets  the  determine  the  relative  strengths  of  the 
head  waves . 
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The  refv’acted  arrival  from  the  layer  Is  eomposec,  of 
head  waves  and  multiply  reflected  wavt?  which  cross  the  .Veyer  at 
least  once  in  the  compresslonal  mode.  All  those  waves  pi-opagate 
with  a  phase  velocity  which  at  large  offsets  approaches  the  com- 
pressional  velocity  in  the  layer.  The  reflected  waves  have  the 
same  shape  as  the  source  function.  The  shape  of  the  head  waves 
la  determined  by  the  time  integral  of  the  source  function.  To 
explain  the  significance  of  the  factor  ("td/CH/Vi  )ywhich  appv.ars 
in  (7)  and  (8)>  we  must  re-examine  the  expression  for  the  direct 
wave  in  (5).  I'he  total  energy  radiated  into  the  system  is 

^  b^/  \l\  } 

'o 

T  is  the  duration  of  the  source  function  (Jjf  )  and,  for  the  form 
of  source  function  we  are  assamlng,  T/i^j  is  a  constant.  The 
total  energy  input  is  directly  proportional  to  "tj .  We  make  the 
energy  input  independent  of  the  dominant  period  by  dividing  the 
source  function  by  the  dimensionless  quantity  (i<l/(^V^i  Th® 
quantity  (i^/(H/'V,  ))*^appears  in  the  numerator  in  (8)  because  the 
3l.L’)e  of  the  head  wave  is  determined  by  the  integral  of  the  source 
function  and  integn’^  ion  introducef*  a  factor  "t^.  Note  that  in 
the  high-frequency  limit,  the  reflected  waves  determine  the  ampli¬ 
tude  and  character  of  the  refracted  arrival. 
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The  complete  expression  for  the  refracted  arrival  is 


When  Is  the  travel  time  of  the  first  reflected  arrival 

which  is  Incident  beyond  the  critical  angle  -  this  is  the  reflected 
wave  P|  P|  .  When  ,  a  shear  head  wave  is  generated  which 

propagates  wj ch  the  phase  velocity  ,  Beyond  the  shear  critical 
angle,  detemined  by  the  travel  time  of  the  shear 

head  wave  - 


t  _  H 
*■"■7 


Inside  the  shear  critical  angle,  is  detemined  by  the  travel 
time  of  the  first  reflected  wave  which  is  incident  beyond  the 
critical  angle  .  Note  that  inside  the  critical  angle 
there  are  no  head  waves  and  if  we  remove  the  restriction  on  f] 
(l.e.,ins:l)  equation  (ll)  gives  the  complete  reflected  arrival 
from  the  embedded  layer. 
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CHARACTERISTICS  OP  THE  MULTIPLY  REFLECTED  WAIVES 

For  high  frequencies  the  multiply  reflected  waives 
deterroine  the  characteristics  of  che  refracted  arrival.  The 
travel-time  curves  for  fh  =  0,1,2  rre  plotted  in  Figure  3*  "The 
difference  in  travel  time  between  the  reflected  wave  and  the 
refracted  arrival  is  plotted  on  the  ordinate.  The  offset  is 
plotted  on  a  logarithmic  scale  to  emphasize  the  fact  that  all 
the  travel -time  curves  for  a  particular  tf\  value  approach  th-^ 
saroo  asymptote  at  large  offsets  -  the  travel-time  difference 
between  any  two  waves  which  have  the  same  m  value  goes  to  zero 
at  larg;^  offsets  regardless  of  the  n  values.  The  travel-time 
difference  between  waves  which  have  different  HI  values  approaches 
a  finite  value  at  large  offsets  which  is  a  multiple  of  T^  .  T^ 
is  defined  in  (9)  v^here  it  is  shown  that  2T^  is  the  time  Interval 
between  head  wave  arrivals . 

This  clustering  of  the  travel-time  curves  leads  us  to 
refer  to  waves  with  the  same  yy}  value  as  a  particular  order.  If 
the  significant  amplitudes  in  each  order  are  confined  to  a  limited 
range  of  w  values  and  T  (source  duration)  ^  ^  ,  the  refracted 
arrival  at  large  offsets  appears  as  a  series  of  events  (Berckhemmer 
and  Oliver,  1955).  Each  event  is  associated  with  a  particular  |7] 
value  and  propagates  with  a  phase  velocity  which  is  very  close  to 
the  compresslonal  velocity  in  the  layer.  The  amplitude  of  each 
e^^rnt  is  obtained  by  summing  the  amplitudes  of  all  rays  which 
have  the  same  fTI  value  as  follows; 
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n 

is  plotted  against  the  normalized  offset  in  Figure  4.  At 
large  offsets,  the  slope  approaches  -3/2.  In  Figure  5f  A^^ 
is  plotted  against  tne  normalized  layer  thickness.  For  thin 
layers  the  slope  approaches  1/2.  A  study  of  many  cases  shows 
that  for  ^/H>-20  and  E/H  -c.  .1, 

(13) 

At  large  offsets  the  head  wave  contributions  decay  like  (p/H)  . 
This  means  that  the  reflected  waves  determine  the  character  of 
the  refracted  arrival  at  large  offsets. 

The  and  depend  on  the  r  ompresslonal  velocity 
ratio,  the  Poisson's  ratios  ( and  0^)  and  the  density  ratio. 

In  Figure  6,  the  7^  and  are  plotted  against  the  compressional 
velocity  ratio  for  a  rigid  layer  (  “  *25)  embedded  in  a  very 

p 

weak  medium  (  0*1  =  .45).  For  large  velocity  contrasts,  the 
'^^(fflil)  approach  and  even  exceed  7^.  'Y^  is  within  10  per  cent 

of  7^  over  a  large  part  of  the  range.  Tj  goes  through  zero 
between  .86  and  .87.  This  fact  indicates  that  in  the  first 
order  (rn  =  1 )  the  amplitudes  of  the  individual  waves  change  sign 
at  least  once.  If  the  amplitudes  of  the  head  wave  and  reflected 


p 

Note  that  the  measurements  of  McDonal,  et  al .  (1958)  indicate 
that  at  a  depth  of  500  ft,  the  Pois.son's  ratio  in  Pierre  Shale 
is  about  .42. 
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wave  contributions  to  the  refracted  arrival  were  detemined  by 

and  alone,  the  reflected  waves  would  be  of  minor  importance 
over  most  of  the  range.  However,  increasing  either  the  offset 
or  the  dominant  frequency  reduces  the  bead  wave  amplitude  rela¬ 
tive  to  the  reflected  wave  No  significance  should  be  attached 
to  the  behavior  of  near  Vj =  1.  This  is  a  consequence  of 
the  fact  that  the  high-frequency  asymptotic  representation  for 
the  head  wave  is  not  valid  when  the  receiver  is  near  the  critical 
cone.  As  Vj /V^  approaches  unity,  the  critical  cone  approaches 
the  receiver  regardless  of  where  the  receiver  is  located.  Note 
that  decreases  relative  to as  the  velocity  contrast 


Inci'eases  and  is  always  less  than  five  per  cent  of  - 

In  Figure  7,  the  andyW^  are  plotted  against 

\  /'Vjj  •  This  quantity  determines  the  Poisson's  ratio 
!  _| 

cr  XJiL_  . 

y±  «! 

TE 

The  higher  oi'^ers  peak  sharply  near  =  3.0(C^  =  .25). 

Near  =  3.0j  is  about  10  per  cent  of  0^  .  As  the 

layer  becomes  more  fluid-like,  the  0^  (rTlil)  decrease  and  become 

very  small  compared  with  Yo  .  This  is  because  the  Yy^  (0)2:1 ) 

depend  on  the  conversion  coefficients  and  R*.,.  As  the  layer 

BS  Sf 

becomes  more  fluid-like,  the  efficiency  of  the  compressional  to 


shear  converslorA  decreases.  Similar  behavior  i'^  exhibited  by  the 
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In  Figure  8  the  '/  and  M  are  plotted  against 

m  m 

Over  the  entire  range  the  (  m^2)  are  small  com¬ 
pared  with  but  as  the  infinite  medium  becomes  more  fluid¬ 
like,  *7,  approaches  Y  and  even  exceeds  it.  The  (m~  2) 

I  0  '  m 

have  much  larger  values  when  the  infinite  medium  is  riuld-llke 

than  when  it  Is  rigid  but  only  ever  exceeds  10  per  cent  of 

A>  •  Note  that  for  large  values  of  all  the  higher 

orders  have  opposite  polarity  to  . 

* 

Equation  (12)  cannot  be  used  when  the  individual  wave¬ 
lets  (corresponding  to  the  different  H  values)  are  wholly  or 
partially  resolved.  The  degree  of  resolution  Increases  a?  the 
offset  decreases,  the  layer  gets  thicker  or  the  dominant 
frequency  Increases.  In  Figure  9,  the  amplitude  of  each  wave 
is  plotted  along  the  ordinate  and  the  difference  between  the 
arrival  time  of  the  reflected  wave  and  the  onset  of  the 
refracted  arrival  along  the  abscissa.  Here  we  consider  only 
multiples  which  do  not  cross  the  layer  in  the  shear  mode  (171  =  0). 
Each  curve  is  drawn  for  a  particular  offset.  The  leftmost  dot 
on  each  curve  gives  the  amplitude  and  time  difference  for  the 
wave  which  crosses  the  layer  twice  in  the  compresslonal  mode 
(n  =  2).  As  w<*  move  along  each  curve  from  left  to  right,  fl 
increases  in  increments  of  two. 

Note  that  the  time  Interval  between  the  arrival  of 
events  cori'e spending  to  different  ITI  values  is  .106.  At  each 
offset,  all  the  significant  arriplitudes  occur  well  within  this 
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interval.  It  follows  that  if  the  first  event  is  resolved  at 
p/H  =  it  ia  resolved  at  all  larger  offsets. 

The  open  dot  on  each  curve  corresponds  to  H  =  8.  As 
the  offset  Increases,  the  amplitude  increases  to  a  maximum  and 
then  decreases.  For  f)  =  8,  the  maximum  amplitude  occurs  near 
f/H  =  20.  At  p/li  =  50,  the  amplitude  for  H  =  8  Is  larger 
than  for  any  o-her  t\  value.  As  the  offset  increases,  the 
amplitude  of  the  later  arrivals  increases  and  the  number  of 
arrivals  within  a  fixed  time  Interval  increases.  When  the  wave¬ 
lets  are  partially  resolved,  this  causes  additional  wiggles  to 
appear  on  the  tall  of  the  first  event  (ftl  =  0 ) . 

As  the  offset  Incrsasep,  the  time  delay  decreases  for 
each  n  value.  The  wavelet  associated  with  each  reflection  moves 
forward  in  time  with  respect  to  the  onset,  and  the  amplitude  of 
the  early  arrivals  decreases.  It  is  this  behavior  vdilch  is 
responsible  for  shingling. 

It  may  seem  odd  that  at  a  given  offset,  the  araplitude 
actually  increases  as  the  number  of  reflections  increases.  ?or 
irn  =  0  the  amplitude  is  determined  by 

A  T  T  R 

At  the  critical  angle  (ft  ),  Rp_  =  1  and  Tq-  =0.  As  0  decreases, 

rr  r,t£ 

Rpp  decreases  rapidly  and  is  always  less  than  unity.  This  means 
that  the  amplitude  reaches  a  maximum  value  near  ft  and  that  the 
position  of  the  maximum  migrates  toward  ft  as  H  Increases. 
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In  Figure  10,  is  plotted  as  a  function  of  Vj /Vg^  -  SU10 

for  several  H  values.  If  the  angle  of  incif'^nce  were  inde¬ 
pendent  of  n  ,  the  amplitude  would  decrease  with  increasing  H 
Actually,  increasing  H  causes  the  angle  of  Incidence  to 
decrease.  This  fact  and  the  proximity  to  the  critical  angle 
of  the  peaks  for  the  higher  H  values  explain  the  increase  in 
amplitude  with  Increasing  ft  .  The  dashed  lines  are  drawn  for 
particular  p /H  values. 
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THE  REFRACTED  ARRIVAL 

The  reflected  waves  and  head  waves  superpose  to  form 
the  refracted  arrival.  The  actual  response  function  is  obtained 
by  (a)  convolving  the  source  function  with  the  sequence  of 
Impulses  which  represent  the  reflected  waves  and  (b)  convolving 
the  integral  of  the  soui*ce  function  with  the  sequence  of 
Impulses  which  represent  the  head  waves  (as  indicated  in 
Equation  (ll)  ) .  The  choice  of  source  function  is  not  completely 
arbitrary.  It  must  be  chosen  so  that  each  poirt  in  an  infinite 
medium  returns  to  its  original  state  a  finite  time  after  the 
arrival  of  the  disturbance  at  that  point .  This  means  that  the 
D.C.  component  of  luch  quantities  as  the  stress  and  strain  must 
be  finite.  The  Laplace  transform  of  the  radial  stress  in  a 
si^erlcally  symmetric  compressional  wave  is 


In  order  for  the  D.C.  compcncnt  to  be  finite. 

Limit  X  -  A .  (I'D 

s-*o  ' 

As  a  consequence  of  (l4),  must  have  at  least  two  axis  crossings. 
We  also  require  that  the  particle  velocity  and  acceleration  be 
continuous.  One  function  which  satisfies  all  these  conditions  is 

)  cos  2^  ,  o-.±^K  J  . 


The  first  factor  determines  the  envelop.  The  duration  of  the 
function  is  T  =  kItj  .  The  condition  2  forces  the  low- 
frequency  behavior  to  satisfy  fl4).  The  source  function  and 
its  integral  are  plotted  in  Figure  11  for  K  =  4.  The  source 
function  is  symmetric  about  't/'t =  2.  The  integral  of  the 
source  function  is  antisymmetric  about  "2.  It  follows 

that  under  certain  conditions,  we  can  distinguish  between 
reflected  and  head  wave  contributions  to  the  refracted  arrival 
on  the  basis  of  symmetry. 

Figures  12-1^  show  how  the  character  of  the  refracted 
arrival  vai'les  with  the  dominant  frequency  at  three  offsets. 

In  each  figure,  the  dimensions  are  iixed.  One  division  along 
the  abscissa  Is  equal  to  the  dominant  period.  As  we  move  down, 
the  dominant  period  Increases,  the  dominant  wavelength  increases, 
and  E/ decreases.  Because  the  dominant  period  changes,  the 
time  scales  cannot  be  compared  directly. 

If  the  refracted  arrival  consisted  of  the  head  wave 
only,  all  traces  wouD.d  have  the  antisymmetric  form  of  the  head 
wave  and  would  t€;rmlnate  at  the  fourth  division.  Figure  12 
shows  three  distinct  events  for  E/,.^  -  1,  5  and  4.  These  events 
correspond  to  rn  =  0  ,  1  and  2.  As  £/>,  decreases,  the  later 

'  '  cj 

events  move  toward  the  first  event,  interfere  with  it  and 
become  lost  in  it. 

Changes  in  E/^^  produce  slgnif'<;ant  changes  in 
character  for  P/H  values  of  3.5  and  10,  but  relatively  minor 
cl.anges  for  P/H  =  100,  At  large  offsets,  the  dominant  frequency 
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is  not  high  enough  to  partially  resolve  the  Individual  waves 
which  contribute  to  the  first  event.  At  small  offsets,  the 
individual  waves  are  separated  by  larger  time  Intervals  and 
tne  change  in  character  is  a  consequence  of  interference. 

As  we  move  down  on  each  figure,  the  dominant  period 
increases.  This  causes  the  individual  waves  to  overlap  to  such 
an  extent  that  only  a  single  event  is  discernible  and  reduces 
the  amplitude  of  the  reflected  waves  relative  to  the  head  waves. 

s 

This  explains  why  the  last  trace  on  each  figure  approaches  the 

antisymmetric  shape  associated  with  the  head  wave. 

Each  trace  is  nomallzed  so  that  the  maximum  amplitude 

plots  at  the  same  value.  The  true  maximum  amplitude  in  the 

refracted  wave  train  is  plotted  in  Figure  15  as  a  function  of 

the  dominant  wavelength  to  thickness  ratio.  Each  curve  is  drawn 

for  a  particular  offset.  As  the  offset  increases,  a  relative 

maximum  develops  which  migrates  toward  the  vertical  line  marked 

"t^  -  .  We  have  shown  that  at  lar^e  offsets  the  waves  arrive 

in  groups  which  are  associated  with  different  m  values  and  are 

=  T.  is  the  condition 
a 

for  constructive  interference  between  the  groups.  The  weak 
minimum  and  maximum  at  p/H  ==  100  are  also  a  consequence  of 
Interference  between  the  groups. 

As  the  dominant  frequency  Increases,  the  individual 


separated  by  the  time  interval  T, 


waves  in  the  group  begin  to  separate  from  one  another,  and  the 
amplitude  varies  in  a  complex  way.  Destructive  Interference 
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between  inulvldual  waves  is  responsible  for  the  deep  mlnlmurr. . 

At  sufficiently  high  frequencies,  the  Iridividual  waves  are 
completely  resolved  and  each  curve  must  approach  an  asymptote 
with  slope  -1/2.  At  very  low  frequencies,  the  head  wave  is 
dominant  and  each  curve  must  a^.  roach  an  asymptote  with  slope 
1/2.  The  increase  In  amplitude  at  low  frequencies  Is  indicative 
of  the  failure  of  the  high-frequency  theory.  If  this  Increase 
in  amplitude  were  continued  to  very  low  frequency,  the  energy 
in  the  refracted  wave  would  exceed  the  total  energy  radiated 
into  the  system  (E.^  ).  There  is  also  another  Indicator  of  the 
failure  of  th«  asymptotic  theory.  We  are  dealing  w'th  a  layer 
which  is  embedded  between  two  half-spaces  which  have  identical 
properties.  In  the  limit  as  the  layer  thickness  goes  to  zero, 
the  am.plitude  of  the  refracted  arrival  must  go  to  zero.  The 
multiply  reflected  waves  exhibit  the  correct  behavior 
(Equation  13 ),  but  the  head  wave  amplitudes  do  not  go  to  zero 
(8).  We  do  not  attach  significance  to  dominant  wavelength  to 
thickness  ratios  greater  than  unity;  consequently,  this  analysis 
is  restricted  to  thick  high  velocity  zones. 

Figures  l6-l8  illustr'ite  the  effect  of  offset.  On  each 
figure,  the  E//^  ratio  is  constant  and  the  dominant  period  does 
not  change  from  trace  to  trace.  Each  trace  is  normalized  so  that 
the  maximum  amplitude  plots  at  the  same  value.  The  nomallzed 
offset  (  p/H)  changes  in  increments  of  0.5.  The  group  of  lines 
which  start  in  the  upper  left  corner  of  Figure  l6  show  how  the 
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individual  peaks  and  troughs  move  forward  with  respect  to  the 
onset  and  d 'crease  in  amplitude  as  the  offset  increases.  This 
behavior  is  known  as  shingling.  The  decrease  in  amplitude 
would  be  considerably  more  pronounced  if  true  amplitudes  had 
been  plotted. 

In  their  model  experiments,  Levin  and  Ingram  (1962) 
observed  that  the  phase  velocity  measured  as  a  function  of 
layer  thickness  tc.xes  on  values  which  exceed  the  compressional 
velocity  in  the  layer.  The  explanation  for  this  effect  lies  in 
the  fact  that  the  phase  velocity  of  the  reflected  waves  exceeds 
the  compressional  velocity  in  the  layer.  The  effect  io  not 
observed  when  the  layer  is  very  thick  because  the  initial  peaks 
and  troughs  are  associated  with  the  true  head  wave  which  is 
observed  uncontaminated  by  reflections  off  the  bottom  of  the 
layer.  For  thin  layers  the  effect  is  not  observed  because  the 
phase  velocity  of  the  reflected  waves  approaches  the  compressional 
velocity  in  the  layer. 

The  secor  i  group  of  lines  in  Figure  18  l-.ow  how  addi¬ 
tional  peaks  and  troughs  appear  on  the  tall  of  the  first  event 
as  the  offset  increases.  These  phases  are  associated  with  the 
later  arriving,  higher  n  values  which  increase  in  amplitude  with 
offset.  The  third  group  of  lines  shows  that  shingling  also  occurs 
in  the  second  event  (fT)  =  l). 

The  shingling  effect  may  be  particularly  pronounced  when 
the  difference  between  the  phase  velocity  of  the  reflected  waves 
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(C)  and  V2  is  large.  The  phase  velocity  is  defined  by 

Cnrh  _  H _  ^  1  . 

Y,  J+nm  SIHG 

'  °wvr 


The  difference  in  phase  velocity  is 

Um-V.  _ 


(15) 


h  JlL  cinQ 

®nm 

Changes  in  the  physical  dimensions  which  reduce  inc  reasc 

the  phase  velocity  difference.  Therefore,  increasing  the  layer 


thickness  or  reducing  the  offset  increases  the  phase  velocity 
difference.  Figures  16-I8  were  computed  for  a  relatively  thin 
layer  (E/H  =  .1).  For  thicker  layers,  individual  extrema  move 
touard.  the  onset  more  rapidly  with  increasing  offset. 

In  Figure  17j  -he  phase  velocity  is  slightly  greater 
for  the  third  peak  than  it  is  for  the  earlier  extrema.  Shingling 
is  exhibited  only  over  a  limited  range  of  offsets.  The  second 
trough  and  third  peak  on  the  first  trace  can  be  correlated  across 
the  entire  record,  but  the  phase  velocity  difference  disappears 
for  p/H  ^  7-  In  Figure  18,  E/^^  =  4  and  there  is  no  shingling. 
These  observations  indicate  chat  for  a  fixed  dominant  period, 
shingling  is  not  observed  beyond  a  certain  maximum  offset.  The 
maximum  offset  decreases  as  tue  dominant  period  Increases  anu 
for  sufficiently  long  dominant  periods,  shingling  does  not  occur. 
When  shingling  is  observed,  it  indicates  that  the  high  velocity 
zone  is  thick.  The  lower  bound  or  E//^^  depends  on  the  depth  to 
the  layer.  For  E/H  =  ,1,  the  existence  of  shingling  Indicates 
that 


Normal  mode  propagation  in  whloh  the  phase  velocity 
exceeds  the  group  velocity  can  produce  shingling  but  not  distance' 
limited  shingling.  This  is  because  t.he  group  velocity  -  phase 
velocity  relationship  for  a  single  mode  is  Independ'^nt  of  distance. 

The  maximu.m  amplitude  and  the  amplitudes  of  individual 
peaks  (P^  )  and  troughs  (T^  )  are  plotted  against  offset  in 
Figure  19 .  Notice  how  the  maximum  amplitude  shifts  to  later 
extrema  as  the  offset  increases.  The  first  extrem'im  may  not  be 
observed  because  of  noise  or  insufficient  dynamic  range.  Nevei - 
thelesSj,  in  certain  cases  it  may  still  be  possible  to  determine 
the  onset  of  the  event.  Figures  12-14  show  that  the  position  of 
the  deepest  trough  is  relatively  insensitive  to  offset  and  to 
E/^jj  when  is  near  unity.  In  each  case,  the  deepest 

trough  occurs  at  a  time  near  2i:^j  after  the  onset.  Recall  that 
is  the  time  at  which  the  source  function  takes  on  its  maxi¬ 
mum  value.  This  suggcists  that  the  time  of  the  source  function 
maximum  might  determine  the  delay  factor  to  be  used  in  determin¬ 
ing  the  true  onset  of  the  refracted  arrival.  When  used  in  this 
way,  the  source  function  should  Include  the  effect  of  reverbera¬ 
tion  in  the  near  surface  in  the  vicinity  of  the  source  and 
receiver. 

In  Figure  20,  the  maximum  amplitude  is  plotted  against 
offset  for  different  values  of  .  Changes  in  the  way  in  which 

the  individual  wavelets  Interfere  accounts  for  the  complicated 
variation  of  the  higher  frequency  curves.  For  E/^  values  of 
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^  anu  I,  tne  amplitude  actual;.y  increases  with  offset  over 

limited  ranges  of  p/H.  For  p/H  ^  8.5,  the  maximum  amplitude 

occurs  for  E/ ==1.  At  laj'ge  offsets  (  P/H  >  20),  each  of  the 

3/2 

curves  approaches  the  (p/H)  decay  required  by  (13). 


CONCLUSIONS 


At  large  offsets  the  refracted  arrival  from  a  layer 
results  from  a  superposition  of  all  waves  which  travel  with 
phase  velocities  equal  to  or  near  the  compressional  velocity 
in  the  layer.  This  requirem.ent  is  satisfied  by  the  head  waves 
generated  by  reflected  waves  which  do  not  cross  the  layer  in 
the  compressional  mode  and  by  multiply  reflected  waves  which 
cross  the  layer  at  least  once  in  the  compressional  mode. 

Specifying  fl  (the  number  of  crossings  in  the  compressional 
mode)  and  fO  (the  number  of  cr-osslngs  in  the  shear  mode)  does 
not  define  a  particular  ray  but  rather  a  family  of  degenerate 
rays  which  arrive  at  exactly  the  same  time.  The  number  of 
I'ays  in  the  family  is 

I 

The  individual  rays  differ  from  one  another  in  the  sequence  in 
which  the  P  and  S  legs  of  the  path  ar*e  traversed.  Because  all 
the  rays  arrive  at  exactly  the  same  time,  only  the  total  ampli¬ 
tude  is  physically  significant .  A  formula  is  derived  which 
gives  the  total  amplitude  in  terms  of  a  series  which  contains 
either  ffl  or  n  +  1  terms  (whichever  is  smaller).  This  formula 
makes  possible  the  rapid  computation  of  highly  degenerate  events. 

The  travel -time  curves  for  all  rays  with  the  same  fH 
value  approach  the  same  asymptote  regardless  of  the  n  values. 

The  travel-time  curves  for  waves  with  different  ffl  values  approach 
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.isymptoteG  which  uiTfer  by  a  multiple  of  .  At  large  offsets 
and  for  thin  layers^,  the  total  amplitude  obtained  by  summing 
all  rays  with  the  same  iTl  value  satisfies 
,  ,  _3/2^  ^  .1/2 

The  head  waves  decay  with  distance  like  (  p/H)~^.  Therefore,  tne 
reflected  waves  predominate  at  large  offsets.  The  reflected  waves 


also  predominate  when  the  dominant  frequency  is  high.  If  the 


duration  of  the  source  function  is  less  than  T^  ,  the  refracted 
arrival  consists  of  a  series  of  events.  Each  event  is  associated 


with  waves  which  cross  the  layer  the  same  number  of  times  in  the 

shear  mode.  At  large  o.ffsets  the  onsets  of  the  events  are 

separated  by  T^  ,  the  phase  velocities  approach  the  compresslonai 

-3/2 

velocity  in  the  layer,  the  amplitudes  decay  like  (p/H)"  an, 
the  relative  amplitudes  are  determined  by  the  .  The  7^ 
depend  only  on  the  density  and  velocity  ratios.  The  higher  orders 
are  shown  to  be  significant  when  the  compresslonai  velocity 
contrast  is  large,  the  Poisson’s  ratio  in  the  layer  is  near  .2b» 
or  the  Poisson's  ratio  in  the  infinite  medium  approaches  ,5. 

The  multiply  reflected  waves  propagate  with  phase 
velocities  which  exceed  the  compresslor.il  velocity  in  the  layer. 
Consequently,  as  the  offset  Increases,  each  wavelet  moves  towarxi 


the  beginning  of  the  refracted  arrival.  In  this  process  the 
early  arrivals  (corresponding  to  the  smaller  H  values)  decrease 
in  amplitude.  Th^  slightly  higher  phase  velocity  and  the  decrease 


-  35  - 


in  amplitude  produce  the  shingling  effect.  The  iatei  arrivals 
(corresponding  to  higher  fl  values)  first  increase  and  then 
decrease  in  amplitude  This  causes  additional  wiggles  to  appear 
on  the  tall  of  each  event  as  the  offset  increases.  The  shingling 
is  distance-limited.  The  maximum  uffcct  at  which  it  is  observed 
decreases  with  increasing  period,  and  for  sufficiently  long  period 
it  is  not  otoervec.  .Shingling  occurs  when  the  high  velocity  zone 
is  thick  compe.red  with  the  dominant  wavelength.  We  know  that  the 
smallest  value  of  (E//ij)  for  which  it  occurs  depends  on  E/H  and 
that  as  E/H  increases,  the  shingling  effect  becomes  more  pro¬ 
nounced. 

For  near  unity,  the  time  of  the  deepest  trough 

is  not  ve^  sensitive  to  changes  either  in  offset  or  in  . 

In  our  case  the  time  between  the  onset  of  the  refracted  wave 
and  the  dee;>'at  trough  is  determined  by  the  time  at  which  the 
source  function  attains  its  maximum  value.  If  thi?'  were  true 
generally,  it  would  provide  a  means  for  making  a  good  estimate 
of  the  onset  time.  When  used  in  this  way,  the  source  function 
must  Include  the  effect  of  reverberation  in  the  near  surface 
in  the  vicinity  of  the  source  and  receiver. 

Beyond  a  certain  offset,  the  maximum  amplitude  of  the 
refracted  arrival  occurs  for  .  This  is  the  condition 

for  constructive  interference  between  the  groups  of  waves 
corresponding  to  different  YT)  values. 
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Our  results  have  been  obtained  by  using  high-frequency 
geometric  ray  theory  and  approximate  the  exact  theory  for 
1. 

Reflections  off  inter'faces  both  above  and  below  the 
high-velocity  zone  can  introduce  energy  into  the  layer  at  and 
near  the  critical  angle  to  produce  secondary  refractions.  A 
reflection  off  a  deeper  interface  can  generate  a  secondary 
refraction  when  the  interface  dips  in  such  a  way  as  to  reduce 
the  phase  velocit,y  in  the  reflected  wave.  It  follows  that  more 
than  one  refracted  arrival  may  be  associated  with  one  high- 
velocity  layer.  When  this  happens,  the  difference  in  travel 
times  between  the  primary  and  secondary^  determines  the  travel 
time  along  the  crit.:.cal  angle  to  an  underlving  or  overlying 


reflector. 
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